Relative phase fluctuations of two coupled one-dimensional condensates 
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We study the relative phase fluctuations of two one-dimensional condensates coupled along their 
whole extension with a local single-atom interaction. The thermal equilibrium is defined by the 
competition between independent longitudinal thermally excited phase fluctuations and the cou- 
pling between the condensates which locally favors identical phase. We compute the relative phase 
fluctuations and their correlation length as a function of the temperature and the strength of the 
coupling. 



I. INTRODUCTION 

Recently, longitudinal phase fluctuations in very elon- 
gated Bose-Einstein condensates have been observed ex- 
perimentally P, 0| ■ Such phase fluctuations are charac- 
teristic of one-dimensional (ID) Bose gases and appear 
in the small interaction regime where p ^Jmpg/Ti, p 
being the linear density of atoms, g the interparticle in- 
teraction between atoms and m their mass. The opposite 
limit, called the Tonks regime Q, where strong correla- 
tions between atoms appear is not investigated in this pa- 
per. For ID Bo se gas es, at temperatures T much smaller 
than T p = Tip^f pg/m/kB, fluctuations of density are su 
pressed and one has a quasi-condensate 0, IE IE 
However fluctuations of phase, given by 
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are still present Q. The logarithmic zero temperature 
term is negligible when using normal experimental pa- 
rameters and phase fluctuations are produced by the 
thermal population of collective modes. 

In this paper we are interested in the case of two elon- 
gated condensates coupled along their whole extension by 
a single-atom interaction which enables local transfer of 
atoms from one condensate to the other. Such a situation 
could be achieved using a Raman or RF coupling between 
different internal states 0- It could also model the case 
of condensates in two very elongated traps coupled by 
a tunnelling effect. The physics of two coupled conden- 
sates, which contains the Josephson oscillations, has been 
studied in a two- mode model in [ IollllllT3 |. In particular 
the many bo dy g round state [Tj and the thermal equi- 
librium state |l2| have been computed. Behind the two- 
mode model the excitation spectrum of two-component 
condensates coupled by a local single-atom coup ling has 
been calculated using the Bogoliubov theory in |l3j . In 
the case of two elongated condensates two effects act in 
opposite directions. Longitudinal phase fluctuations in 
each condensate tend to smear out the relative phase be- 
tween the two condensates, while the coupling between 
the condensates energetically favors the case of identical 
local relative phase. The goal of this paper is to deter- 
mine the relative phase of the two condensates at thermal 
equilibrium as a function of the strength of the coupling. 



II. FORMALISM 
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FIG. 1: Situation studied in this article. Two elongated con- 
densates are coupled by a interaction which enables local trans- 
fer of atoms from one condensate to the other. 

We are interested in pure ID condensates where the 
temperature, the interaction energy and the coupling 
strength 7 are all much smaller than the transverse con- 
finement energy. Thus the Hamiltonian is written 



H = d 



2m 



Q2 Q2 



^(#W+^(^.W]}, (i) 



where ipa,b are the boson annihilation operators for the 
condensates labelled a and b, U(z) is the trapping po- 
tential and /i is the chemical potential. Assuming that 
the size of the transverse ground state a± = yj21i/muji_ 
is much larger than the s-wavc scattering length a, the 
effective coupling constant is simply g = ^l. |f- . 

Following calculations made for ID condensates 0, Q 
we expand the field operators in terms of their density p 
and phase 9 as 



1>a,b(*) = e i9 - b{z) y/paA*)- (2) 

The Hcrmitian density and phase operators obey 
[Pi(z),6j(z')} = i8(z-z')S hJ 0. 

As we are interested in temperatures small enough to 
be in the quasi-condensate regime, density fluctuations 
are small and we write 



Pa,b( z ) = Po( z ) + 8pa,b{ z ), 



(3) 



where 



Po 



■C 1 and po satisfies the Gross-Pitaevskii 
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equation modified by taking p — > // + 7: 



and for the antisymmetric family it becomes 



2m 



A + U{z)-{p + 1 ) + g oPo 



po~=0. (4) 



We also assume that the phase difference between the 
condensates at a given position is small 



|A0(z)| = \9 a (z) - e b (z)\ « 1. 



(5) 



The Heisenberg evolution equations for p a b and 9 a ,b are 
developed to first order in Sp a ,bi ^&a,b and AO and we 
obtain 



hd t O a . b = 
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The first terms on the right hand side are identical to 
those for a single condensate and the second terms cou- 
ple the two condensates. We perform a canonical trans- 
formation to the bosonic operators 
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D„..b — n , + 1\/P0"a,b, 
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which evolve according to 
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where we have introduced the operators 
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Such an evolution is the same as the one given by the 
standard Bogoliubov theory and we recover indeed the 
same result as that of . As the matrix in Eq. © is in- 
variant by exchange of a and b, eigenvectors may be split 
in two families: the symmetric eigenvectors invariant by 
exchange of a and b and the antisymmetric eigenvectors 
which arc multiplied by -1 by exchange of a and b. The 
eigenvalue equations are thus reduced to two 2x2 ma- 
trix equations. For the symmetric family the eigenvalue 
equation becomes 
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As for the standard Bogoliubov theory the Hamiltonian is 
then written, up to a real factor, as a sum of independent 
bosonic excitations 



H 2 = e skb\ k b sk + Y tnkb ] nk b nk 



(13) 



k k 

and the B operator is written 



B a .b = J2( bskUsk + b lk<k) ± J2(bnku nk + bl k v* nk )(lA) 

k k 

where the sums are done only on the eigenvectors nor- 
malized to / dz(\u k \ 2 - \v k \ 2 ) = 1/2. 

We are interested in the correlation function of the 
phase difference AO which is written, after commuting 
the B operators to normal order, 

(AO(z)A0(z')) = (: A6(z)A6(z') :) + fcfl. (15) 

The second term accounts for the phase fluctuations in a 
coherent state with linear density po for a and b. We are 
not interested in this term, and thus we will consider only 
the normal ordered expectation value. If we expand this 
in terms of the b operators and consider thermal equilib- 
rium where no correlations between different excitations 
exist, we obtain 

{:A6{z)A9{z>):) = 1 £ {(SL^X/nVnV + fnk'fnk* 
Po l 



~ v *nkfnk V nkf nk '} 



(16) 



where the prime means that we evaluate the function 
at z' and f~ k — u nk — v nk . As expected only the anti- 
symmetric modes contribute because we are interested in 
phase difference. This expression gives the relative phase 
fluctuations once the modified Bogoliubov spectrum of 
Ea. l|12(l has been calculated. In the following we will 
give explicit results in the case of an homogeneous gas. 



III. RESULTS FOR HOMOGENEOUS 
CONDENSATES 

We now consider an homogeneous gas with periodic 
boundary conditions in a box of size L. The potential U 
then vanishes and the Gross-Pitaevskii equation gives 

p = gp - 7. (17) 

The Bogoliubov function can be looked for in the form 



u sk = (2L) 2 exp(ikz)U sk 
v sk = (2L)~5 exp(ikz)V sk 



(18) 
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where |C/ s fc| 2 — \V s k\ 2 = 1 and similarly for the antisym- 
metric modes. The Bogoliubov eigenvalue equation for 
the symmetric branch then reduces to the standard Bo- 
goliubov equation 



h 2 k 2 

2m 
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whose spectrum and eigenvectors are well known. For 
the antisymmetric case the eigenvalue equation becomes 



gpo 



gpo 




which is simply the same as the symmetric case, with the 
kinetic energy shifted by 27. Thus the eigenvalues and 
eigenvector components are 



Vik 
2m 



'- + 2j+2gp 



U n k + V n k 
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This two-branch spectrum was already obtained in a 
more general case in 01 • In the case where 7 3> gpo, 
these excitations are almost purely particles with V n k -C 
U n k for any k and their energy is simply h 2 k 2 /2m + 27 as 
expected for a particle in the state (\a) — \b))/\/2 and of 
momentum k. In the opposite case where 7 <C gpo, three 
zones can be identified. For k <C 2^/rfvy/Ti we obtain col- 
lective excitations with V ~ U and with energy 2^/jgpo- 
For 2 y /m-f/h, <C k <C 2y/mgpo/fo we still have collective 
excitations with V ~ U but their energy is given by the 
normal Bogoliubov dispersion law hk\J gpo/m. Finally 
for k 2^/mgpo/h excitations are just particles with 
energy H 2 k 2 /2m. 

Using the plane wave expansion i|18|) and the normal- 
ization condition U 2 k — V 2 k = 1 the correlation function 
(|16|1 of the relative phase fluctuation is written 



will be smaller than the phase fluctuations of a single 
condensate which implies 



AQ 2 a V m gPo ln ( Ly/mgp 



hp 



(23) 



The whole theory is valid only for large density so that 
y/mgpo/ {hpo) <C 1 and in the experiments accessible un- 
til now the size of the condensate is not large enough to 
produce noticeable phase fluctuations at zero tempera- 
ture. 

Phase fluctuations are thus due to thermal excitation 
of the collective modes and we will give a simplified ex- 
pression by making several approximations. First we will 
approximate the Bose factor by 



»k = 



k B T 

f-nk 



(24) 



This is justified as this expression deviates in a signif- 
icant way from the Bose occupation factor only when 
nk becomes smaller than 1, ie when e n k > k B T, and 
the contribution to phase fluctuations of those modes is 
small even with the previous expression which overesti- 
mates their population. Let us now consider separately 
the case where 7 3> gpo and the case 7 <C gpo- 

If 7 > gpo, then (U nk - V nk ) 2 ~ 1 for all k and e k ^ 
h 2 k 2 /2m + 27. This gives, approximating the discrete 
sum by an integral, 
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As we consider only temperatures k B T <C Tipo y/ gp jm so 
that we have quasi-condensate, these phase fluctuations 
are always very small. 

Let us now consider the case where 7 <C gpo- The 
modes with \k\ 3> ko = ^/mgpo/h give a negligible contri- 
bution to the phase fluctuations. Indeed for those terms 
{U n k — V n k) 2 — 1 and tk — h 2 k 2 /2m so that their contri- 
bution to the phase fluctuation is 



(: A0(z)A6(z') 



2p L 



{( U nk ~ Vnk f (2n nk + 1) - l} cos[fc(2 - z'%22) 



where n n k = l/(e €nk / kBT — 1) is the occupation num- 
ber for the state with energy e n k- Using the expression 
(|21|l this correlation function can be computed numeri- 
cally. In the following we analytically compute the phase 
fluctuations using some approximations. 

The terms which do not involve rife correspond to the 
zero temperature contribution. As the function V 2 k — 
UnkVnk is always smaller than the corresponding function 
for a single condensate, the relative phase fluctuations 



mk B T 
irh 2 p 



' dk h 



k B T^Jm 
nftpo^/gpo 



(27) 



which is always small in the regime of quasi-condensates. 
Thus only the modes \k\ <C 2^gpom/h are considered for 
which 



(Unk - VnkY 



^Jh 2 k 2 /m + 47 



(28) 



and the correlation function then becomes 



A0(z)A0{z') 



2k B T 



POTT Jo 



dk 



h 2 k 2 



4 7 



cos[fc( 



0129) 



The integral can actually be extended to infinity as higher 
k values give negligible contributions and we find 



A0(z)A0(z') 



k B T 

2p h 



m 
— exp 

7 



(30) 



Note that this expression is the same as Eq.(26), which 
was not expected a priori. This formula, which give the 
amplitude of the relative phase fluctuations as well as 
their correlation length 1/7 is the main result of the 
paper. It agrees well with the numerical calculation of 
Ea. (|22fl as shown in FigEl Phase fluctuations are small 
only if 



fc B T«p ?H/-. 

m 



(31) 



Note that as we assumed small relative phase difference, 
this is also the limit of validity of our calculation. The 
phase diagram of Fig [31 summarizes the previous results. 
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FIG. 3: Phase diagram for the fluctuations of the relative 
phase between the two condensates. Only temperatures much 
smaller than fipoWgpo/ (kB\frn) are relevant as for larger tem- 
peratures one does not have a quasi- condensate anymore. For 
temperatures larger than hpo/(kBLy / 7n), each condensate has 
longitudinal phase fluctuations. Below the curve, which cor- 
responds to Eq. 13 ID . the coupling between the condensates is 
large enough to suppress local relative phase fluctuations be- 
tween the two condensates. Above this curve, there are local 
relative phase fluctuations between the two condensates. 




oscillation frequency is |ToL ITTl IT3 | 
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FIG. 2: Correlation function of the relative phase fluctuations. 
The solid line is the numerical calculation of Eq. \2°A) with 
7 — gPo/W, T = hpo-y/j I '(2y/mkB) and L — 100?i/ y/mgpo. 
The dotted line is the analytical expression Eq. \HU\) which only 
differs from the numerical expression at small separations. 



IV. DYNAMICAL INTERPRETATION 

The condition (|31[) to have small relative phase fluctua- 
tions has a dynamical interpretation which is shown very 
qualitatively below. In a two-mode model of the Joseph- 
son coupling between two condensates of N atoms it has 
been shown that if 7 <C Ndp/dN, then the Josephson 



fip V~fgpo 



(32) 



On the other hand a single elongated condensate will 
experience phase fluctuations and the phase at a given 
position will evolve in time. If the change of the phase 
during a Josephson oscillation time is small, then the 
Josephson coupling will ensure that the relative phase 
between the two condensates remains zero: there will be 
no relative phase fluctuations. However if the change of 
the phase during a Josephson oscillation time is large, 
then the Josephson coupling will not have time to adjust 
the phase of one condensate with respect to the other: 
there will be relative phase fluctuations of the two con- 
densates. We thus have to compute the change of the 
local phase of a single condensate 



((»(iM)-«(o)) : 



(33) 



with the average corresponding to the thermal equilib- 
rium and the coupling between the two condensates be- 
ing ignored. This calculation could be done rigorously 
by developing the operator 9 on the collective excitation 
bosonic operators bk and b k . In the following we present a 
simpler argument that gives the same order of magnitude. 
We first estimate the amplitude Ak of the phase modula- 
tion of wave vector k. The energy of this phase modula- 
tion is just the kinetic energy N\Ak\ 2 h 2 k 2 / 4m which in a 
classical field theory at thermal equilibrium corresponds 
to an energy of k B T/2 and thus 



\A, 



2mk B T 1 
h 2 k 2 ~pvL' 



(34) 
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This is indeed the contribution of the mode k to phase 
fluctuations as computed in Eq.|J2SJ if 7 = 0. According 
to the Bogoliubov spectrum and because only modes with 
|fc| <C 2^/gpom/h contribute, the mode k evolves with the 
frequency 



(35) 



The evolution of the phase after a Josephson oscillation 
time tj ~ Ti/ \Jl9Po is then written, after averaging over 
the independent phases of the phase modulations, 



({8{tj)-6{0)f) 



E 

k 



\A k \ 2 (l-cos(Lu k tj)) 

2mk B T f°° 1 - cos{^/tJgp^k/^/m) 
h po2ir 
ksTy/rh 



,1k 



(36) 



Small relative phase fluctuations of the two condensates 
occurs when this quantity is small and we recover the 
condition (pTTf) . 



V. DISCUSSION 



about 2tt from one end of the condensate to the other as 
soon as T >T (j> = Ti 2 p /{mLk B ) = 1.8 nK. However the 
local relative phase between the two condensates stays 
much smaller than 1 if T <C hpo^/j/{kB^rn) = 180 nK. 
The calculations made here for homogeneous condensates 
could be used to describe a trapped inhomogeneous gas 
via a local density approximation similar to that used in 
|l5j as long as both the healing length lh = h/y/mgpo 
and the correlation length of the phase fluctuations are 
much smaller than the extension of the condensate. In 
the above example, 1^ = 0.6 pra and l c = 2 /im are indeed 
much smaller than L. 

To measure experimentally the relative phase fluctua- 
tions and their correlation length, one should perform an 
interference experiment. In the case where the two states 
are internal states, an intense 7r/2 pulse has to be applied. 
Measurement of the local density of atoms in the state 
I a) and \b) then gives access to the local relative phase 
of the two condensates. In the case where \a) and \b) are 
confined in the wells of a double well potential, the in- 
terference measurement is performed via a fast release of 
the confining potential followed by a time of flight long 
enough for the two clouds to overlap. Indeed, the total 
intensity presents fringes in the direction orthogonal to 
z|l6f and, at a given z, the position of the central fringe 
gives the value of the local relative phase. 



In conclusion we have shown that as long as the tem- 
perature is small enough to fulfill Eq. I|31|) , although there 
might exist large phase fluctuations along each conden- 
sate, the local relative phase of the two condensates 
stays small. In the opposite case there are large fluc- 
tuations of the relative phase whose correlation length 
is l c = h/2^/mj. As an example let us consider the 
case of two Rubidium condensates of 10 4 atoms elon- 
gated over L = 200 pym, confined transversely with an 
oscillation frequency loj_/2tt = 1 kHz and coupled using 
7 = 50 Hz. The phase of each condensate changes by 
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